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ABSTRACT 



The mass shifts of experimentally well-known baryons due to meson-baryon self- 
energy loops are calculated, and their impact on the observed splitting of the baryon 
spectrum is studied. Configuration-mixed wave functions adapted from a 'relativized' 
model are used with the ^Pq model to provide predictions for the strength and 
analytical momentum dependence of the strong vertices. Intermediate states include 
all the lightest pseudoscalar and vector mesons and corresponding baryons required 
to provide a complete set of spin-flavor symmetry related baryon-meson states. The 
sum over intermediate-state baryons is extended to include the second {N — 3) band 
of negative-parity excited states, to provide the most complete calculation of its kind 
to date. 

It is found that with reduced-strength one-gluon-exchange interactions between 
the quarks, roughly half of the splitting between the nucleon and Delta ground states 
arises from loop effects. The effects of such loops on the spectrum of negative-parity 
excited states arc also studied, and it is found that the resulting splittings arc sensitive 
to configuration mixing caused by the residual interactions. With the extensive 
set of intermediate baryon-meson states used, a reasonable correspondence is found 
between model masses and the bare masses required to fit the masses of the states 
extracted from data analyses. 
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CHAPTER 1 
INTRODUCTION 

The field of nuclear physics spans a very wide range of phenomena and, despite a 
long history, is still riddled with a large number of unanswered questions. Within it, 
the field of hadron physics finds itself in a 'bridge' position between the visible world 
of atoms where electrons and nuclei rule, and the realm of high energy physics where 
quarks, gluons, and other such 'invisible' elementary forms of matter are the focus of 
attention. Hadrons are ultimately what our world is made of, i.e. the fundamental 
pieces of matter bound together to become the neutrons and protons (two examples 
of hadrons) , then further combined to form the nucleus of each atom charted on the 
periodic table of elements. 

Studying hadrons therefore means trying to understand how and why particles 
like quarks and gluons are combined to produce the experimentally observed spectra 
of objects labeled 'mesons' (predominantly quark + anti-quark + gluons states, such 
as pions) and 'baryons' (predominantly three-quark + gluons states such, as protons 
and neutrons) as opposed to any other possible combination. Additionally, while 
an electron can be removed from an atom and isolated to study how forces act 
upon it, the 'strong force' that holds quarks together is so strong that we have 
not, and according to our current understanding will not, see a single quark in 
isolation. Matter can be disintegrated down to its quark constituents in highly 
energetic accelerators but the length of time quarks remain individual particles (i.e. 
not bound to others) is incredibly short, and ultimately they end up bound to each 
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other in hadrons. Therefore much about the strong force and how it binds quarks 
into hadrons remains poorly understood. 

Not only is the study of the strong interactions pushing the limit of current 
experimental facilities, but theoretical tools that are known to work well with 
electromagnetic or high energy phenomena have very limited use in the study of 
quarks and gluons, at least at the energy scale where they form bound states such 
as baryons and mesons. In addition, the theory which is now accepted to be correct 
for strong interaction physics is so complex that the equations governing it can only 
be solved exactly for a very limited number of cases, leaving sizeable gaps in our 
understanding. 

Theoretical work is therefore ongoing in two main areas. One is to broaden the 
range of applicability of exact solutions of the theory of strong interactions via novel 
computing techniques and algorithms. The other is to develop and/or improve models 
which, although not exact, give a good overall picture of many manifestations of the 
strong force and open windows to the physics behind many experimentally observed 
phenomena. Within this context, the broad goal of the work presented in this thesis 
is to remove a level of approximation in one already successful baryon spectroscopy 
model, i.e. a model which describes and predicts the number and properties of baryon 
states that should be 'seen' experimentally. Doing so should improve the ability of 
this model to explain some of the intricacies of the baryon spectrum, while increasing 
its predictive power and helping resolve some of the discrepancies between predictions 
and observations. 

1.1 Quantum Chromo dynamics 

It is now accepted that the building blocks of matter are the quarks, leptons, and 
gauge bosons. Quarks come in six different species or 'flavors', and are known as 
the up down 'd', strange 's', charm 'c', top 'f , and bottom '6' quarks. The six 
flavors of leptons are the electron 'e""', muon '/x~', and tau 'r~' and their associated 
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neutrinos Vei^iJ,, and u^-. Quarks and leptons form a group of particles called fermions. 
Each has half-integer spin and an associated anti-particle, for example the positron 
'e"*"' or the anti-top quark which are identical in mass but different in charge and 
color (another quark quantum number discussed more below) to their particles. The 
gauge bosons are integer-spin particles and are the photon '7', Z, W^, gluon 'g\ and 
the proposed graviton. 

There are four fundamental interactions, each with one or more particles associ- 
ated with it that carry the force. Gravity, the weakest of these interactions, is thought 
to be mediated by massless, spin-two gravitons. The weak interaction, whose force 
is carried by spin-one, self-interacting Z and bosons, is really a component of 
the electro-weak force and is responsible for radioactive beta decay processes. The 
other part of this force manifests itself via electromagnetic interactions, which act 
between electrically charged objects and are mediated by spin-one, electrically neutral 
photons that do not interact among themselves. Finally, the strong interaction, 
which acts between color-charged objects, is mediated by self-interacting vector 
bosons called gluons and is responsible for nuclear binding and the interactions of the 
constituents of the nuclei. The quarks, which come in three colors and have fractional 
electric charge, are known to have strong, electromagnetic, weak, and gravitational 
interactions. The leptons, such as the electron, are subject to all forces except the 
strong one, as they do not carry a color charge, while the neutrinos have neither 
strong nor electromagnetic interactions. 

The strong force is 'weak' at very short distances, creating what is called 
asymptotic freedom (quarks appear to behave like free particles) , but grows infinitely 
strong at 'large' distances, resulting in the phenomenon known as confinement (a 
quark cannot be isolated like an electron or a proton) . The strong interactions provide 
the 'glue' to hold the quarks together to make hadrons, the strongly interacting 
particles, the mesons and baryons. Although the number of quarks within hadrons 
is not defined by quantum field theories, several models have baryons composed 
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of three 'constituent' quarks, which have half-integer spin and obey Fermi-Dirac 
statistics. Mesons on the other hand are described as integer spin quark- antiquark 
pairs that obey Bose-Einstein statistics. To be observable, particles such as hadrons 
are required by confinement to be colorless objects (color singlets). Baryons must 
therefore have one red, one blue, and one yellow quark (as RBY — white) and mesons 
are allowed to be a colorless combination of RR, BB, and YY quark- antiquark pairs. 
There is growing evidence for the existence of other quark and gluon states such as 
glueballs (pure gluon states) and hybrids {eg. qqqg), as well as hypothesized states 
such as diquonia {qqqq), dibaryons {qqqqqq), and others. 

The equations describing the electromagnetic interactions were formulated by 
Maxwell and form the basis of Quantum Electrodynamics (QED). Correct theories 
for weak and strong interactions came later. It is now widely accepted that a field 
theory equivalent to QED exists for strongly interacting particles and is known as 
Quantum Chromodynamics (QCD), with gluons mediating the forces between colored 
quarks, which are analogous to photons. The QCD Lagrangian has the form 

C = ~F,,aFr + ^{ip- m)*, (1.1) 

where = d^-igA'^^t'' is the covariant derivative, F^" = d^'A'^-d''A^+g J2 fabcA'^A"^ 
is the field tensor, and are the gluon fields (a = 1 to 8). The last term in the 
field tensor definition indicates that, unlike the photons of QED, gluons interact with 
each other, giving QCD its non-abelian behavior (for more information on QCD see 
for example Refs. [1-3]). Unfortunately, unlike QED, there is as yet no obviously 
successful way to go from the QCD Lagrangian to a complete understanding of the 
large number of observed hadrons and their properties. Lattice QCD, a field theory 
that replaces space-time with a lattice of discretely spaced points (colored sources 
-quarks- at the junctions and color electric fiux lines -mediated by gluons- as links 
between them), is making visible progress towards that goal by using numerical 
techniques to solve otherwise unreachable problems, but calculations beyond masses 
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and static properties of the ground states and the hghtest negative parity excited 
baryon states are still in the future. Other methods, such as expansions based on 
the large Nc (number of colors) limit of QCD, and effective field theories (theories 
that replace part of unknown interactions by physically sound approximations), are 
limited in their scope and are currently unable to provide the global yet detailed 
understanding needed for a description of all aspects of strong interaction physics. 

1.2 Non-perturbative QCD and QCD-based Models 

It is relatively easy to treat electroweak interactions via perturbation theory (PT), 
but strong interactions of hadrons involve dealing with QCD, where one cannot 
expect much of PT in a situation that is fundamentally not one of weak coupling, as 
is explained below. 

The usual way to treat local interactions is through PT, i.e. by expanding 
various quantities in powers of the coupling constant. In QED, it is useful to 
define an effective coupling constant a{Q^), which gives the momentum transfer (Q^) 
dependence of the renormalized vertex function. This function receives contributions 
from vacuum polarization graphs which describe the electron loop corrections to the 
photon propagator. The result is an effective coupling a{Q^) which increases with 
Q^. In QCD however, this situation is complicated by the gluon self- interactions 
(gluons carry color charge). The effective quark-gluon vertex can be summed over 
all orders of the renormalized coupling (see figure 1.1) and has the form 

12n 

"^^^'^ = i33-2Nf)lniQyAlany ^^'^^ 
where Nf is the number of flavors, and Aqcd is the momentum scale at which 
becomes strong as is decreased. 

Prom Eq. 1.2 it can be seen that at large Q^, or short distances (relative to 
the Aqcd scale and for A^^ < 16), agiQ"^) — > (a property known as asymptotic 
freedom) so that hard processes (processes calculable using perturbative QCD as 
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Figure 1.1. Perturbative expansion of the gluon-gluon interactions defining the 
strong coupling constant as{Q^). 

the momentum transfer is large enough to produce a small value of as) such as deep 
inelastic scattering, some weak decays of heavy flavors, and some experimentally seen 
but theoretically highly suppressed strong decays of heavy quarkonia, can successfully 
be treated perturbatively. On the other hand, at small momentum transfers, or 
large distances, as{Q'^) grows quite large and PT becomes invalid. Unfortunately, 
this is the region relevant to strong decays of hadrons composed of light quarks, 
electromagnetic transitions, and the weak decays of hadrons containing light flavors, 
like hyperons (baryons containing one or more strange quarks) or kaons (mesons 
containing one strange quark or anti-quark). Experimental measurements yield a 
value of Aqcd ~ 200 MeV. The strong interactions become strong at distances larger 
than ~ 1/ Kk, 1 fm, which is roughly the size of the light hadrons. 

Within QCD, there are two main phenomena that are essentially non-perturbative 
and therefore cannot be obtained even by summing entire perturbative series. The 
first is confinement, and the second is the dynamical breaking of chiral symmetry. 
Confinement is related to the interaction energy, which increases with distance in 
contrast to the Coulomb energy. The breaking of chiral symmetry gives a dynamical 
mass of several hundred MeV to light quarks, which are known to have current-quark 
masses of only a few McV at large momentum scales. Both phenomena arc connected 
to non-zero vacuum expectation values, which would vanish at any order of PT: the 
so-called quark (< qq >) and gluon (< gg >) condensates. 

There are two main alternatives to PT based on field theory. One is the QCD 
duality sum rules approach, which uses a short distance expansion of products of 
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operators. This method allows certain quantities to be calculated in terms of the 
condensates < qq >, < gg >, and perturbatively calculable coefficients. These 
quantities can also be calculated in terms of hadron properties such as masses, 
widths, and branching ratios, hence the duality. The other is lattice QCD, which 
attempts to completely calculate hadronic properties from first principles (the QCD 
Lagrangian on a lattice of space-time points based on the work of Wilson [4] , see 
Ref. [5] for an introduction to lattice QCD) based on the physical idea that the 
long-distance properties of QCD are the most important (confinement and dynamical 
breaking of chiral symmetry), and that short distance properties can be reached by 
extrapolation. These techniques permit strong coupling calculations based on the 
assumption that the unrenormalized coupling constant is very large, which use Monte 
Carlo simulations based on Feynman path-integrals. 

Recently a covariant approach to the description of hadron structure has been 
developed, based on the Bethe-Salpeter equation and the Schwinger- Dyson method 
of solving field theory. This approach has been widely developed for the description 
of meson, but so far restricted to the descriptions of ground-state baryons. 

Although these methods are increasingly useful, they remain technically very 
complex and somewhat limited in their applicability. Additionally, using these 
techniques, it is often difficult to obtain physical insight (especially on the lattice) 
about phenomena such as decay mechanisms or confinement (although confinement 
is a natural consequence of the lattice, i.e. the area of the Wilson loop gives an energy 
oc r so produces a linearly rising potential, it is not a proof of its existence). This 
highlights the continued need for other methods which arc inspired by QCD but not 
necessarily derived from it. The inability to calculate with QCD in the low regime 
has made it necessary to use phenomenological models of hadron structure based on 
expectations of the low energy behavior of QCD. The quark potential model and 
other dynamical models were created to fulfill this purpose. 
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Within the framework of the quark potential model, the spectra of mesons 
and baryons, as well as their strong, weak, and electromagnetic decays have been 
successfully calculated. This model allows for direct calculation of relevant matrix 
elements for each definite decay and provides transparent direct links to experimental 
data. Its simplicity implies a lack of theoretical foundation in QCD, but despite that 
fact its past and present empirical successes are impressive. 

1.3 Corrections to the Quark Model 

In QCD there are qqq{qq) configurations possible in baryons, and these must 
have an effect on the constituent quark model, similar to the effect of unquenching 
lattice QCD calculations. These effects can be modeled by allowing baryons to 
include baryon-meson (BM) intermediate states, which lead to baryon self energies 
and mixings of baryons of the same quantum numbers. A calculation of these effects 
requires a model of baryon-baryon-meson (BB'M) vertices and their momentum 
dependence. It is also necessary to have a model of the spectrum and structure 
of baryon states, including states not seen in analyses of experimental data, in order 
to provide wave functions for calculating the vertices, and to know the thresholds 
associated with intermediate states containing missing baryons. 

The goal of the present work is to self-consistently calculate such effects for a set 
of experimentally well known baryon states. The method and results arc presented as 
follows. In Chapter 2 the work of several authors who have made contribution to this 
field is reviewed, and important elements are extracted and related to the present 
work. The nature and extent of the present research is then discussed in more 
detail, highlighting the improvements needed to be made to this type of calculation. 
Chapters 3 and 4 present an overview of the different methods adopted for use in this 
research. Results are then presented in detail in Chapter 5, and finally Chapter 6 
offers conclusions and outlook for future extensions of this work. 
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CHAPTER 2 
BARYON SELF ENERGIES 

Baryon self energies due to BM intermediate states and BM decay widths can 
be found from the real and imaginary parts of loop diagrams (see Figure 2.1). The 
size of such self energies can be expected to be comparable to baryon widths. For 
this reason, they cannot be ignored when comparing the predictions of any quark 
model with the results of analyses of experiments. Since the splittings between states 
which result from differences in self energies can be expected to be comparable to 
those that arise from the residual interactions between the quarks, a self-consistent 
calculation of the spectrum needs to adjust the residual interactions, and with them 
the wave functions of the states used to calculate the BB'M vertices, to account for 
these additional splittings. 

Earlier studies have brought these facts to the attention of the nuclear physics 
community, each highhghting different aspects of the problem. What seems to be 
missing is a consistent and complete calculation of the effects of the self energies 
on the baryon spectrum. By looking back at the existing literature, lessons can be 
learned on how to accomplish such a task as thoroughly as possible. 

2.1 Existing Work 

Previous calculations of the self energies of ground state and negative-parity 
excited baryons use baryon-meson intermediate states consisting of ground states. 
The work of Zenczykowksi [6] takes the point of view that the 'residual' interquark 
interactions are unimportant, and that hadronic loop effects dominate the observed 
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Figure 2.1. Contributions from a baryon-meson loop to a) tiic self energy of baryon 
B and b) the mixing of baryons with same quantum numbers. The dashed arc 
represents a virtual meson. 

splitting and mixing pattern of the ground and first (negative-parity) excited states 
of baryon states in the octet and decuplet SU(3) / representations. 

This work claims that at least two thirds of the observed splittings in these states 
can be attributed to such effects, as can the mixing angles between states due to 
these effects. In a simplifying limit, a formula relating the E — A and A — mass 
differences, derived using onc-gluon exchange by de Rujula, Georgi, and Glashow [7], 
can be attributed to the effects of hadronic loops. 

This calculation uses only spatial ground state intermediate baryons, but unlike 
those of some other authors, considers a complete set of accessible SU(3) / interme- 
diate states. This means that for and A baryons the intermediate states with 
pseudoscalar mesons were A"7r, Att, Nt], At], Nt]', At]', SAT, AAT, S*Ar, and those 
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with vector mesons were Np, Ap, Ncu, Acu, EK*, AK*, and 'E*K*. This calculation 
also considered the self energies of strange baryons, where a different (and larger) 
set of intermediate states is possible. A dispersion integral relates the shift in the 
squared mass of a baryon state to properties of the intermediate states and a spectral 
function m'g,mM), which depends on the nature of the baryons B (initial) and 
B' (intermediate) involved in the loop diagram, 

ms-{ms) ^l^w, /(,)— -2 ^ds. (2.1) 

% tnr i> 

Here the sum runs over all (open and closed) decay channels i = B'M with 
\/s^j. — niB' + ruM, and the weights wf give the spin and flavor [SU(6)] overlaps 
between the B and B'M states. For B and B' baryons restricted to ground states, 
and without SU(3)/ breaking in these ground-state baryon wave functions, these 
weights have the important property that 

E^f = 48, (2.2) 

i 

for all baryons B, as long as the sum runs over all intermediate states i allowed by 
the quantum numbers. This means that in this symmetry limit the mass shift due to 
these effects is the same for all of the ground state baryons, and so no mass splittings 
are generated by loop effects, as might be expected. This observation is critical, as 
it makes clear that without the inclusion of at least this set of intermediate states, 
calculations of these effects do not start from the symmetry limit and so cannot be 
expected to give physically meaningful results. 

Away from this limit the nature of the spectral functions p becomes important. 
Their values were calculated using the ^Pq model with universal radii for the mesons 
and baryons, taking into account the spin, flavor, and spatial structure of the baryons 
and mesons involved. 

The author concludes that, with reasonable hadron radii, about two thirds of 
the splittings and mixings in ground and negative-parity excited state baryons must 
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arise from hadronic loop effects, and that this allows the use of a significantly smaller 
coupling constant ct^ in the quark residual interactions which may explain the rest of 
the splittings. It is crucial that a 'complete' set of SU(6)-related B'M intermediate 

states is employed. 

The work of Blask, Huber, and Metsch [8] is similar to that of Zenczkowski, except 
that vector mesons are not taken into account in the intermediate states, and the 
point-like mesons are coupled to baryon states using an elementary-meson emission 
(EME) model with the meson-quark coupling fixed from the ttNN, KAN and rjNN 
coupling constants derived from analyses of experimental data. The recoil energy 
of the intermediate baryon is neglected, but intermediate state thresholds are given 
by using the physical masses of the intermediate hadrons. Baryon masses including 
these loop effects are found by diagonalising an effective Hamiltonian which includes 
a term for the internal dynamics of the baryon and meson, and an energy-dependent 
term 

5H,{E) = Hc ^^l_^ Hc, (2.3) 

where He describes the coupling between baryons and mesons and H is the full 
Hamiltonian. 

Similar conclusions about the importance of such effects are made, although this 
calculation suffers from a point- like treatment of the mesons (leading to overestimated 
widths for decaying baryon states) and some problems in the resulting spectrum 
which hkely arise from the restricted set of intermediate states. This calculation also 
hints at a possible cancellation between the spin-orbit effects due to the one-gluon- 
exchange residual and splittings which arise from the inclusion of B'M intermediate 
states. 

Brack and Bhaduri [9] calculate self energies of the nucleon and A ground states 
only, using only pions as intermediate mesons, but do not restrict the intermediate 
N and A baryon states to spatial ground states. They find that the difference in the 
self energies of the nucleon and A ground states converges to within 5 MeV of the 
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large N result when a set of intermediate baryons up to and including the N — 3 
band (second band of negative-parity orbital excitations) states is employed. They 
find that, in their model, the difference in the pionic self energies of the odd-parity 
excited states and the ground state converges too slowly to make definite statements. 

Part of this trouble with convergence may be due to their model of the BB'M 
amplitudes, which simply attaches a pion to the quarks with a (nonrelativistic) 
pseudoscalar coupling, with an additional axial form factor 

F.(k^) = l/(l + kVA^), (2.4) 

with — 1275 MeV, corresponding to the mass of the ai meson. Since their loop 
amplitudes involve elementary intermediate pions, they include a factor of 1/uJk-i 
where ujk = \J^^ + rn^ is the pion energy, from the normalization of the wave function 
of the intermediate pion. This factor is not present in the pion center of mass wave 
function in nonrelativistic models which treat it as a composite particle. Although 
the presence of this factor has the effect of further suppression of high-momentum 
contributions to the integral over the loop momentum, the net result is that it is still 
likely that the effective pion-nucleon vertex in this model is too soft. In subsequent 
models and the present work a more rapid decrease of the vertex amplitudes with 
is shown to produce better results for the mass shifts, and can be attributed to an 
effective size for the operator which creates a constituent quark-antiquark pair (see 
Geiger and Isgur [10]). 

The intermediate states arc described by simple unmixed harmonic oscillator wave 
functions. The excitation spectrum of the intermediate states is taken to be either 
harmonic oscillator plus zero-range (contact) spin-spin potential, or the Isgur-Karl 
potential which is modified by anharmonicities in the spin-independent potential, 
which gives a more realistic spectrum for the energy of the intermediate states. They 
show that, at least for the miclcon-A splitting, the details of this spectrum are 
unimportant. It can be expected that they will become very important, however, if 
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this calculation was extended to the self energies of excited states, as these depend 
crucially on the positions of the thresholds due to the opening of various channels 
for excited states to decay to excited states. 

The important conclusions of this work are: convergence of the nucleon-A mass 
difference in the sum over excited intermediate states can be demonstrated; it is 
misleading to include only baryon ground states as intermediate states, as inclusion of 
excited states reduces the difference in the nucleon and A self energies substantially; 
their final results depend sensitively on the chosen (axial) radius of the nucleon, as 
expected, and changing the gluonic hyperfine splitting changes the difference in the 
self energies of the nucleon and A; and if the gluonic hyperfine splittings are too large 
(> 250 MeV) it is impossible to fit the observed A — N splitting. It is also noted that 
one-pion exchange can, with some adjusted parameters (a reduced strength coupling 
to the quarks), be made to simulate these effects. Poor convergence was found in the 
calculation of the self energies of the negative-parity excited states, an issue which 
will be resolved in the present work. 

The work of Horacsek, Iwamura and Nogami [11] would appear to partly con- 
tradict that of Brack and Bhaduri, with A — = 20 MeV from the inclusion of 
baryon-pion intermediate states. However, the approximation of each intermediate 
state quark moving in a single central potential (shell model) is used, so that inter- 
mediate excited baryon states are described as individual excited quark substates. 
This means that the intermediate states are far from a basis of hadrons, and so 
this calculation ignores what we know about the spectrum of confined hadrons in 
the intermediate state, and the resulting thresholds. Both of these calculations are 
incomplete because they do not include contributions from mesons other than pions. 

Silvestre-Brac and Gignoux [12] examine the self energies of only the lowest 
lying negative-parity excited states, and focus on total spin 1/2 and 3/2 spin-orbit 
doublets in the A^, A, A, and S fiavor sectors. They correctly use a complete set of 
SU(6)-related intermediate states, but as in Zenczykowksi's work, these are restricted 
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to spatial ground states. No configuration mixing is allowed between states due to 
interquark Hamiltonian (Hq), so spin-orbit partners in these negative-parity excited 
states are degenerate. The calculation also uses only one radius for all baryons and 
all mesons, i.e. uses the simplifying assumption of SU(3)/ symmetry in the wave 
functions. Bare masses are solved for self-consistently, i.e. are free parameters. The 
decay thresholds are found using physical masses for the intermediate virtual hadrons, 
which is equivalent to adopting dressed states in the propagators of the intermediate 
hadrons, and this is of course crucial to obtaining a correct description of widths. 
The calculation also uses a cut-off factor in the integrals over the loop momentum. 
The authors justify this by comparison to calculations which use elementary pions 
in the intermediate state and so have a further suppression factor of the inverse of 
the pion energy, 1/cUk, and from a lack of information about strong vertices at large 
relative momenta of the final-state hadrons. 

Their conclusions are that hadronic loops are important ingredients in the 
understanding of spin-orbit splittings, with a satisfactory description of the order 
and magnitude of the spin-orbit splittings of negative-parity excited baryons resulting 
from their calculation. However, this latter conclusion seems premature given that 
it has been shown by Brack and Bhaduri [9] and Gcigcr and Isgur [10] that the 
restriction of the intermediate state baryons to ground states results in self energies 
which have not converged. 

The calculation of Pujiwara [13] uses antisymmetrized (3g)(gg) cluster-model 
wave functions composed of simple harmonic oscillator wave functions and the 
plane-wave relative motion to describe the baryon-meson intermediate states. The 
decay operator employed is unlike those in other calculations, as pairs are created by 
an interaction between a quark and a quark- antiquark pair creation vertex which is 
consistent with the residual interactions between the quarks in the hadrons (see also 
Ackleh, Barnes and Swanson [14]). In particular, it contains the contact, tensor, and 
spin-orbit interactions arising from one-gluon-exchange between the quarks. The self 
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energies of ground states and lowest lying negative-parity excited states of N, A, 
A and E baryons are calculated using intermediate states restricted to ground state 
pseudoscalar and vector mesons, and ground state octet and decuplet baryons. The 
nonrelativistic approximation is made in the energy denominators in order to allow 
analytic treatment of the loop integrals involved in the evaluation of the self energies. 

Rather than adopt the momentum dependence of the vertex amplitudes which 
arises from the structure of the hadron states, the simplifying assumption of a 
universal dependence on the relative momentum k of the intermediate hadron 
pair is adopted. As in other calculations, this dependence is modified from the 
exp{—k'^/6a,^) dependence given by a nonrelativistic evaluation of the decay ampli- 
tude in the presence of recoil, where a is the harmonic oscillator size parameter, in 
order to further suppress high-A; contributions to the loop integrals. In this case it 
is simply given the value exp(— A;^/3a^). Flavor-symmetry breaking is ignored in the 
pair-creation interaction for simplicity. 

The results show that it may be possible to arrange a cancellation between 
spin-orbit splittings arising from the interactions between the quarks and from loop 
effects, and to describe the mixings and decay widths of these states in the same 
model. Notable exceptions are the flavor singlet (lowest lying) negative-parity A 
states A(1405) and A(1520) which are about 100 MeV too heavy, as in simple 
three-quark models. 

As mentioned previously, conclusions made in the models described above about 
spin-orbit forces in negative parity excited baryons are likely to be premature, given 
the information provided about convergence by Brack and Bhaduri [9]. It is shown 
in the present work that the inclusion of negative-parity excited baryons in the 
intermediate states and configuration mixing in their wave functions are crucial to 
the accurate calculation of mass shifts of these states. 

From the above work it is clear that a self-consistent and successful model of 
baryon self energies must employ a complete set of spin-flavor symmetry related 
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B'M intermediate states, and at the same time must include excited baryon states 
up to at least the N — 3 band in order for the sum over intermediate states to have 
converged. This requires a detailed and universal model which is capable of relating 
the baryon spectrum and the decay amphtudes of a wide variety of baryon states 
to a wide variety of baryon-meson final states in an efficient way. It is also clear 
that it will be necessary to modify the usual momentum dependence of the decay 
amplitudes calculated in this model to take into account the size of the constituent 
quark-pair creation vertex. 

In addition, the size of these loop effects requires that the interactions between 
the quarks required to fit the observed spectrum be changed by the presence of these 
loop effects. It is inconsistent to not then also change the wave functions used to 
calculate the vertex amplitudes and to examine the effect of these changes on the 
self-energies. Brack and Bhaduri [9] have shown that the A-nucleon splitting may 
not be sensitive to such details, but from the sensitivity to the structure of the 
interquark Hamiltonian used to describe the hadron states observed in many of these 
calculations, it can be expected that this will be an important effect in the calculation 
of the self energies of the negative-parity excited baryons. 



Based on all the lessons learned above, our goal is then to calculate the energy- 
dependent self energy of baryon B given by 



where Mbb'm is the analytical strong decay matrix element of initial baryon B 
decaying into two hadrons, baryon B' and meson M, as calculated with the ^Pq pa-ir 
creation model. The integral is taken over the relative momentum k between baryon 
B' and meson M. 



2.2 Current Work 




(2.5) 
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Eq. 2.5 comes from first-order, time-ordered perturbation tlieory, and is not a 
fully relativistic, frame-independent equation. Here it is evaluated in the rest frame 
of the decaying baryon B. It has both real and imaginary parts but the real part 
only is extracted by taking the principal part of Eq. 2.5, in effect performing the 
integration everywhere along the real line except over a small symmetric interval 
centered over the pole location. 

We calculate the self energies of Eq. 2.5 for the ground states Nucleon and A, and 
lightest negative-parity excited states. As detailed further in appendix A, baryons 
have color, spin, flavor, and spatial wave functions obtained from the different ways 
three colored, flavored, spin- 1/2 quarks can be combined with the relative angular 
momenta of the three-body system's two relative coordinates (see flgure A.l) p and 
A. With no orbital angular momentum, i.e. L = Ip -|- ^ = and Up — n\ — 0, 
only two states are possible: = 0^ ® {S = ^ oi S = ^} ^ = ■, f"*"} 
where the total angular momentum J = L -|- S and parity P = (—1)'''"'"'^. The 
hghtest of these states is the Nucleon, with — and flavor wave function 
either uud (proton) or udd (neutron). The other combination gives us the A with 

= I"*" and flavor wave function uuu (A"*""*"), uud (A"*"), udd (A°), or ddd (A~). 
The next lowest-lying states come from the addition of one unit of orbital angular 
momentum {Ip — 1 oi Ix — 1, and Up — n\ — G). The states form the flrst 
band [N — 2{np + n\) -\- Ip -\- l\ — 1] of negative-parity excited states produced 
via — 1~ ® {S = \ oi S = ^} ^ — {\ , I , I } to give the seven states 
2 iV| ; A| ; 2 ; A| ; which are constructed based on the overall symmetry 
of the combined wave functions. 

There are two main ingredients needed to complete such a calculation. The flrst 
is a model of the spectrum and structure of baryon states. The second is a model of 
baryon-baryon-meson vertices and their momentum dependence. 

The model of the spectrum must include not only states seen in analyses of 
experimental data, but also states classifled as 'missing', in order to provide wave 
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Figure 2.2. Momentum structure of the strong decay vertex A4 BB'M{k) in the initial 
baryon B center-of-momentum frame. 

functions for calculating the vertices, and to know the thresholds associated with 
intermediate states containing missing baryons. As mentioned above, the splittings 
between states which result from differences in self energies can be expected to be 
comparable to those that arise from the residual interactions between the quarks. 
A complete calculation of the spectrum therefore needs to adjust the residual 
interactions, and with them the wave functions of the states used to calculate 
the baryon-baryon-meson vertices, to account for these additional splittings. The 
relativized quark potential model used in this work is presented in more detail 
in Chapter 3, including the changes required by the presence of the additional 
baryon-meson intermediate states. 

The formalism used to model the strong decay vertices is presented in more detail 
in Chapter 4. It provides an analytical form of the momentum dependence of each 
vertex M.BB'Mik) as a function of the relative momentum k between the intermediate 
baryon B' and meson M in the center-of-momentum frame of the initial baryon B 
(see Figure 2.2). Based on the lessons learned from earlier work, this calculation 
includes all allowed combinations of intermediate states B'M from the sets 

M e {tt, K, rj, rj', p, co, K*}, B' e {N, A, A, E}. (2.6) 
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including all excitations of the baryons states up to and including the N — 3 band 
states. Excited mesons states have been omitted at this time as their higher mass 
and additional angular momentum highly suppresses decays to these states. Studies 
including more massive initial states should, however, consider including such states. 

Based on the range of intermediate states included, it is worth noting that for 
each A'"^*) initial state studied, the sum of Eq. 2.5 {J2b'm) includes a total of 591 
intermediate baryon-meson states. Similarly, a total of 378 intermediate states are 
included for each A^*) initial state. 

As the self energies due to a given intermediate state depend crucially on the 
masses adopted for the intermediate hadrons, these are taken to be the physical 
masses, where known, and model masses [15] otherwise. The 'bare' mass required to 
reproduce the known physical mass of any initial baryon state B considered is found 
by solving the self-consistent (highly non-linear) equation 

E + ^b{E) = Mb (2.7) 

for the 'bare' baryon mass E'^. Therefore the integration of Eq. 2.5 needs to be 

performed over a range of bare masses Eb in order for the final result {E'^) to be 
extracted from the intersection of the left side of Eq. 2.7 with the right side when 
Mb — Mphysicai- 
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CHAPTER 3 
THE QUARK POTENTIAL MODEL 

The nonrelativistic constituent quark model (NRQM) owes its origin to many 
authors but the model of Isgur and Karl and collaborators [16-19] has been quali- 
tatively successful in both meson and baryon sectors, despite its lack of theoretical 
foundation in QCD. In an effort to correct some of the flaws of the NRQM, the 
Isgur-Karl model was later 'relativized' by Godfrey and Isgur [20] (for mesons) and 
Capstick and Isgur [21] (for baryons). This latter version has been used by Capstick 
and Roberts [15, 22-23] in extensive calculations of strong decay amphtudes. A 
modified version of this model is used in the present work to obtain the masses and 
wave functions of known and 'missing' baryon states. It is therefore appropriate 
to give an overview of the main components and discuss the value of some of the 
parameters used here. 

The choice of dynamical degrees of freedom used to represent a baryon depends 
on momentum transfer. At low Q^, they can be taken to be constituent quarks, which 
are valence quarks with effective masses of about 220 MeV for u and d (~330 MeV in 
the NRQM), and about 420 MeV for the s quark (~550 MeV in the NRQM). In this 
model the gluon fields affect the quark dynamics by creating a confining potential 
in which the quarks move. At short distances, a perturbative one-gluon exchange 
between quarks is assumed to provide the spin-dependent potential. 

3.1 The Hamiltonian 

The Hamiltonian used [21] for the baryon system has the form 
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i/=i^O + K>ge + Konf. (3.1) 

where Hq is the relativistic kinetic energy term 

Ho = j:{pl + mlY/', (3.2) 

i=l 

Voge is the one-gluon-exchange potential, and Konf consists of a string potential and 
its associated spin-orbit term arising from the Thomas precession. 
The one-gluon exchange potential has the form 

i<j 

with the color induced interactions being 

y-oge yCoulomb , y-hyperfine , yspin-orbit(cm) /g ^\ 



where 



yhyperfine y contact _j_ y tensor (3 5) 



The Coulomb term is spin-independent and proportional to l/vij (where rjj is the 
relative position of the (ij) quark pair), the spin-orbit and hyperfine interactions are 
color-magnetic in nature, and the hyperfine interaction consists of a Fermi contact 
term oc S^{rij) and a tensor piece. The terms of the one-gluon exchange potential 
can be found from the Brcit-Fermi reduction of the one-gluon exchange T-matrix 
element oc u{p', s')7^h(p, s) where u is approximated as the Dirac four-spinor of a 
free particle. 

The confining potential is composed of two parts 

Vstring ^b^Vij, (3.6) 

and 

T/ _ T^spin-orbit{Tp) fo 7\ 

^^spin-orbitrs) — / . ^ii ■ W-'j 

i<j 

The string part of Konf is the adiabatic potential corresponding to the energy of 

the minimum-length configuration of the Y-shaped string linking the quarks. The 
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spin-orbit term includes the Thomas precession effects of the full spin-independent 
potential. For ease of calculation, V^tring is approximated by a sum of a constant, an 
effective two-body piece, and a three-body piece 

V^tring = Cqqq + + ^3b, (3-8) 

where Cqqq is an overall energy shift which arises from the vacuum modifications due 
to the presence of colored fields in the baryon, / = 0.5493 is chosen to minimize 
the size of the expectation value of Vsb in the harmonic oscillator ground state 
of the baryon system, and h is the meson string tension. The two-body part of 
Ktring is calculated directly during the diagonalization of the Hamiltonian, and Vsb 
is computed perturbatively. 

The potentials have been modified from their nonrelativistic limit {p/m — > 0) 
by several effects. For example, since constituent quarks are not point-like, the 
interquark coordinate is smeared out over mass-dependent distances. The smearing 
is brought about by convoluting the potentials with a function 

A,(r,,) = 4e-'^^-^. (3.9) 
7r2 

where the (Tij are chosen to smear the interquark coordinate over distances of 
approximately 0.22 fm for light quarks, and 0{1/Mq) for heavy quarks Q. A second 
modification allows the potentials to be momentum dependent by introducing factors 
which replace quark mass terms by energy dependent ones such as 

miiTLj 

^ ip!j + mn'/'ip!j + m])y^ ^ ^^^^ 

where Pij is the magnitude of the momentum of either quark in the ij center-of-mass 
frame. These terms are included in the potentials in the form of factors such as 
where the e^s are free parameters designed to allow the rough description 
of the momentum dependence of each potential. 

23 



3.2 The Parameters 

For completeness, we include the final value of some of the relativized quark 
potential model parameters used in this work. As will be explained in more detail in 
Chapter 5, the values used in Ref. [21] and subsequent work were not adequate here as 
they were meant to model a spectrum without taking into consideration the existence 
of qqq + qq configurations. Therefore our work requires slightly different parameters. 
Another reason for reducing the value of (Tq, which corresponds to the inverse of a 
quark 'size' and is one of the smearing parameters used to define (Tij in Eq. 3.9 (the 
other parameter being s), is that it brings the electromagnetic form factor of the 
quark required to fit nucleon electromagnetic form factors in relativistic (light-cone 
based) models more in line with this strong size. Studies have been done with various 
values of some of the parameters to understand their effects before selecting the final 
values. Some of the 'intermediate' results will be presented in Chapter 5 to illustrate 
this process. More information about the different potentials, the origin and use of 
the parameters listed here can be found in Ref. [21] and references within, since their 
description is beyond the scope of this work. 

It is important to note that all spin-orbit effects have been removed from the 
Hamiltonian used to obtain the wave functions used in this work. Studies on how 
best to introduce spin-orbit effects are in progress but are secondary to the main goal 
of this work, and so are not presented at this time. 
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Table 3.1. The parameters of the relativized quark potential model. 



Parameter 


This work 


Ref. [21] 




^{mu + md) (MeV) 


220 


Same 


Light quark mass 


rris (MeV) 


419 


Same 


Strange quark mass 


b (GeV^) 


0.15 


Same 


String tension 


2 ~^ ^cont 


1 - 0.168 


Same 


Relativistic factor 


2 ~^ ^tens 


1 - 0.168 


Same 


u 


1 + ^Coul 


1 

2 


Same 


u 


^critical 


0.550 


0.60 


as{Q' = 0) 


ao(GeV) 


0.833 


1.80 


Relativistic smearing 


s 


1.55 


Same 
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CHAPTER 4 



MODELING STRONG DECAYS 

4.1 Introduction 

One very important element of our calculation is a model of the momentum 
dependence of each vertex in Figure 2.1 (a). If that diagram were to be cut in half, one 
could see that each half represents a decay B — > B'M. We can therefore use a decay 
model to obtain the structure of the vertex and hence its momentum dependence. 
The two most suitable classes of models for this work are briefly described below 
before going into the details of the specific model used here (for a recent review of 
these and other strong decay models see Ref. [24]). 

The first class, known as elementary-meson-emission (EME) models, has baryons 
treated as objects with a quark structure while mesons arc treated as elementary, 
point-like objects emitted from a quark during the decay. Each decay transition is 
described in terms of a coupling constant. This implies many parameters, although 
several coupling constants can be approximately related via SU(2) or SU(3) flavor 
symmetry. This class of models lends itself well to relativistic treatment, which is 
often desirable for light mesons such as the pion. Unfortunately since mesons are 
modeled as point-like objects, treatment of excited mesons is restricted since radial 
excitations imply an extended spatial wave function which is not modeled. 

The other class of models, referred to broadly as pair creation models, treat all 
hadrons as composite objects. The decay of a hadron coincides with the creation of 
a quark-antiquark pair somewhere in the hadronic medium. The created antiquark 
then combines with a quark of the original hadron to create a daughter meson 
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while the created quark becomes part of the other daughter hadron. These models 
describe hadron emission in a unified way, often involving only one free parameter 
(the pair-creation strength 7) , and allow for the treatment of all excited baryons and 
mesons within the same framework. In contrast to EME models, pair creation models 
are non-relativistic and therefore approximations are made. They are nonetheless 
more realistic, simple, and have been successfully applied to the study of a broad 
range of strong decays for both mesons and baryons. 

There are several types of pair creation models, such as the ^Pq and ^5"! models 
(named after the quantum numbers of the created pair), and the flux-tube and string 
breaking models, where the location of the created pair is restricted to an area inside 
the chromoelectric fiux tube (the tube of gauge field that is shown by lattice QCD to 
form between two colored sources) or along the string axis. We describe our choice 
in some detail in the next section. 

4.2 The Model 

Due to its simplicity and past successful applications to the strong decays of 
hadrons, the ^Pq model, popularized by Le Yaouanc et al. [25] has been selected to 
be used in this research. It has been widely applied to baryon decays [15] [22] [23], 
meson decays, and even generalized to the decay of states composed of n-quarks [26]. 

Within this model, the strong decay can be seen as a process where a quark- 
antiquark pair is created from the QCD vacuum with quantum numbers J^*^ = 0"^"*". 
As shown below, in the '^^'^^Lj notation, this corresponds to ^Pq, hence the name of 
the model. The pair can be created anywhere in space, but wave function overlaps 
will naturally strongly suppress creation very far from the initial hadron. The created 
pair is added to the initial system, giving rise to two new non-interacting final state 
hadrons. To be observed, these new hadrons must be color singlets. Additionally, 
the qq pair must be neutral with respect to the additive quantum numbers, meaning 
that it must also be a fiavor singlet, and have zero total angular momentum. Because 
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Figure 4.1. OZI-allowed (a) and suppressed (b) processes B B'M, in a quark 
pair creation model. 

the quark and antiquark have opposite intrinsic parity, parity conservation further 
dictates the pair be in a relative p-wave (i.e. L = 1) so that its total spin must be 
one {S — 1) to combine to the required J — 0. 

It is important to note that only Okubo-Zweig-Iizuka (OZI) allowed decays are 
considered, as illustrated in Figure 4.1(a). A process is said to be OZl-forbiddcn, or 
suppressed, [see Figure 4.1(b)] if the quark of the created quark pair does not combine 
with quarks in the initial hadron but instead the created quark and antiquark form 
a separate meson. 

4.2.1 The Operator 

The starting point in modeling the B B'M transitions of baryons in the '^Pq 
model is the form of the operator T responsible for the decay. Within this model, 
the operator does not result from a detailed Hamiltonian that would come from the 
QCD Lagrangian, as the complexity would be overwhelming. Instead it is entirely 
phenomenological and is defined only for the decay process under consideration. It 
has the following form 
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X ^(l,m;l,-m|0,0) xT, 3^r™(Pi -P.) ^1(pO 4(Pi)' (4-1) 

m 

where Cjj and Fij are the color and flavor wave functions of the created pair, 
both assumed to be singlet, Xij is the spin triplet wave function of the pair, and 
3^1 (Pi ~ Pj) = |Pi ~ Pjl ^i(Pi ~ Pj) is the solid harmonic indicating that the pair is 
in a relative p-wave (L = 1). Note that the threshold behavior resulting from this 
I Pi — Pjl factor is as seen experimentally. Here 6|(Pi) and d}j{Pj) are the creation 
operators for a quark and an antiquark with momenta Pi and pj respectively. The 
exponential has been introduced to give the vertex a spatial extent by creating the 
quark- antiquark pair over a smeared region, instead of at a point as is the case in 
the usual version of the ^Pq model. The addition of this form factor 'softens' the 
vertices and suppresses the self energy contributions from intermediate states where 
the hadrons have high relative momentum. 

There are only two phenomenological parameters in this model. The first one is 
7, the coupling strength, which we fit to the experimentally well known A — > Nn 
decay, and the second one is /, which is set to give a reasonable quark-pair-creation 
vertex size of around 0.35 fm (the same as that used in Geiger and Isgur [10] and 
Silvestrc-Brac and Gignoux [12]). 

Gonsider an initial observable system A (a baryon composed of three quarks) 
decaying into two observable, non-interacting hadrons; baryon B and meson C. One 
quark from A will merge with the created antiquark to form meson C, and the 
remaining two 'initial' quarks will merge with the created quark to form baryon B. 
The notation used is illustrated in Figure 4.2. Note that in this version of the ^Pq 
model, quarks 1 and 2 are considered spectators as they do not participate in the 
decay. 
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A:Sa^Jp, + 1/2; 

Ja = Sa + La„ 



C : = 1/2 + 1/2; 




Figure 4.2. Schematic diagram of the decay B B'M in the ^Pq modeL 
The angular momentum notation is shown. The decay proceeds through 
5(123) ^ 12(44)3 ^ 5'(124)M(43). 

4.2.2 Wave Functions Considerations 

For the transition A — > EC, we are interested in evaluating the following 
transition amplitude 

M^{BC\T\A), (4.2) 

where \A) denotes the wave function of the initial baryon A, and \BC) the wave 
function of the final baryon-meson pair. The initial system is assumed to be 
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in a static state (made up of a quantum superposition of harmonic oscillator 
substates) , multiplied by a plane wave with relative state momentum for center- 
of-mass motion, where a denotes the quantum numbers needed to describe the basis 
states. The total wave function for the initial state 1^4) expressed in momentum 
representation is expanded in terms of basis states 

= 5(p„ - K„)*^ = ^d^^f'^^ (4.3) 

a 

The coefficients d^a ai'c obtained by diagonalization of the Hamiltonian H in 
the basis of the , taken here to be harmonic-oscillator basis states. The series is 
truncated to = 2{np + nx) + lp + lx — 6 for the positive-parity states and = 7 for 
negative-parity states, giving of the order of 100 substates for each J^. Note that the 
level of the expansion is related to the sum of the powers of the two coordinates 
p and A (see Appendix A) in the associated Lagucrrc polynomials. A higher power 
means a shorter length scale, therefore the maximum N was chosen to 'resolve' the 
shortest range interaction in the Hamiltonian, or equivalently allow the variational 
calculation of the energy to converge. The expansion coefficients are such that 
the total wave function \['^'^" is antisymmetric, despite the fact that the basis states 
are taken to be antisymmetric only in the first two quarks. 

The wave functions for the final baryon and meson are given in a similar fashion 

= 6{p, - K,)*^ = d^K"^'- (4-4) 
h 

^^.K. ^ ^^p^ _ ^^^^c ^ d^^^'^^ (4.5) 

c 

By combining ^^'^t and we obtain a wave function \BC] JbMhKbJcMcKc) 

which describes the hadrons in a plane wave with their angular momenta decoupled. 
For ease in further treatment of angular momenta, we couple J^-I-Jc = Jbc, and change 
the variables and Pc to K = -|- Pc and k = | (p^ — Pc) . Then K represents the 
total momentum of the BC system and is conserved through the term S{pa — K^) of 
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equation 4.3, and k is the relative momentum between B and C. Instead of using a 
plane wave |k), we change to a spherical wave llmk) via 

(p6Pe|Ko/mA;o) =5{K- Kq) —^5{k - ko) . (4.6) 

Finally, the relative momentum / is coupled to Jbc to give the total angular momentum 
Ja, so that the form of the final state wave function becomes 

\BC, JbJc, Jbcl'i JaMa, Kq/Co) 

= / dKhdK^ Yl {JbMbJcMc\JbcMbc) 

Mbc,m,Mb,Mc 

x{J,,M,Jm\JM{K,K,\Kolmko)\BC,JbM,K,J,M,K,). (4.7) 
Baryon states are written as 

^^Ca<I>J2^X- (4.8) 

where Ca, 4>, ip, and x color, flavor, spatial, and spin wave functions 

respectively. The baryon wave functions used in this calculation were produced using 
a relativized model [21] with variable-strength spin-dependent (one-gluon exchange) 
contact, tensor, and spin-orbit interactions between the quarks. More details about 
the baryon wave functions can be found in Appendix A. 

4.2.3 Transition Amplitude 

From equation 4.1, the transition amplitude is not Galilean invariant since it 
contains a factor 5p, where p is evaluated in a definite frame. The results therefore 
depend on the chosen frame of reference. A good choice of frame is the one where 
the decaying baryon A is at rest, so we set = 0. Momentum conservation yields 
a factor 5(Ko) in the amplitude, and we now rewrite equation 4.2 as 

{BC\T\A) = S{K,)Ma^bc. (4.9) 

Incorporating equations 4.3, 4.4, and 4.5, we obtain the expression 

Ma^bc = E di^*d';:*dtMA^Bc{a, b, c) (4.10) 

a,6,c 
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The color, flavor, spin, and spatial degrees of freedom can be separated via invariance 
techniques. The resulting amplitude then reduces to products of sums over internal 
summation variables of 6—j and 9—j coefficients from angular momentum recoupling, 
and flavor, and spatial matrix elements that can be calculated independently, which 
then can be combined to give a total decay matrix element. The final form of 
MA-^Bc{a, b, c) is 

"^V<3 Jp,Sa,Sb 

( Sa Lp "Sa 1 f Lp Sp Lp) ( Sb Lp "Sft 1 f Lp Sp Jp 1 

14 Ja LJXI Sa SaiXh Jb Lb } \ ^ Sb Sb f 





'Jp 


1 
2 


Sb 




' Sb 


h 


Jb 




1 

2 


1 
2 


Sc 




Sc 


Lc 


Jc 


Sbc 




1 


Sbc_ 




_Sbc 


Lbc 


Jbc 



xE^'{'£ si t}{1' j! t}^(4,i^c,i^6c,^,4,i^,A;o), (4.11) 

where the factor of 6 comes from the redefinition, for the created pair, P = pj + 
and p = |(Pj — Pi) so that the spherical harmonic found in the operator (eq. 4.1) can 
be rewritten as — 873^1 (—2p) = 673^1 (p). The overall factor of | is the color matrix 
element, JF is the flavor overlap, and e is the spatial matrix element. 

Further explanation of notation and derivations of some components have been 
gathered in appendix B for the interested reader. 
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CHAPTER 5 
RESULTS 

Graphical and tabulated results for the self energies of the ground state Nucleon, 
A and non-strange L — 1 negative-parity states are presented in the following 
sections. First, evidence of the convergence of our results is shown, indicating that a 
minimum number of intermediate states have been included in order to obtain stable 
and reliable results. The same graphs also show the effects of the various decay 
thresholds, and their effect on the sums of self energies. Next, tables displaying the 
impact of changes in the Hamiltonian and the associated baryon model wave functions 
are presented. Finally the qqq spectrum obtained from the modified Hamiltonian is 
graphically compared with the spectrum of bare energies obtained from fitting the 
sum of the bare energies and self energies to the physical masses. This illustrates 
that it is possible, in a self-consistent calculation, to describe the observed masses 
with a combination of splittings induced by interquark forces and differences in the 
self energies. 

5.1 Convergence and Thresholds 

One important result coming out of this calculation is the phenomenon of 
convergence. As pointed out before, the number and type of intermediate states 
included in this type of calculation can dramatically change the final results. In the 
figures that follow, this concept and the associated consequences will be illustrated 
for the states studied. 

Figures 5.1 through 5.6 show the sclf-cncrgy contributions to the masses of several 
baryons, from the sum of intermediate B'M states for B' including progressively 
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higher harmonic oscillator bands, and M the complete set of pseudo-scalar and vector 
mesons {i.e. tt, K, rj, rj', p, cu, K*). These figures show the results obtained from using 
wave functions created from a Hamiltonian that includes both the contact and the 
tensor part of the hyperfine interaction. In each figure, each line represents the sum 
E + 'Eb{E) with J2b{E) obtained with a set of baryons B' in the harmonic oscillator 
bands indicated in the subscript [T,{E)n=o, E(£')jv<i, E(£')jv<2, and I!,{E)n<s\. For 
each of these sums of intermediate states, the corresponding bare mass E^ can be 
extracted by reading the value of the energy E corresponding to the intersection 
between the curve for the sum of self energies E + I1b{E) and the horizontal line 
representing the physical mass. This process is in effect solving 

E + Eb{E) = Mb, (5.1) 

for E — E^ when Mb — -^physical with progressively larger sum over bands of baryon 
intermediate states A^ = 0, A^<1,A^<2, and < 3. 

As will be seen in some of the figures and tables that follow, occasionally more 
than one solution is possible for Eq. 5.1 due to oscillations caused by the presence of 
B'M decay thresholds (the energies at which the decays B B'M become physical, 
i.e. energetically possible). In these few cases a range of values is presented unless it 
is clear that one solution is favored. More details are presented in the next section. 

5.1.1 Nucleon and Delta Ground States 

Figure 5.1 and Table 5.3(b) show the self energy contributions to the mass of the 
Nucleon from a progressively larger sum of intermediate baryon-meson states. The 
first point to notice is the sizeable difference between the 'bare mass' E'^ when only 
baryons in the N — band are included, i.e. Eg — 1.85 GeV, and when other baryon 
states are included, Eq — 2.36 to 2.50 GeV. These bare masses are not observables, 
but the mass sphtting between the nucleon and other states is, therefore any large 
variation in the bare mass of the nucleon can change its relationship with other states 
in the spectrum. 
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Next note that this 'bare mass' difference is not 'built' of equal contributions from 
each set of intermediate states but comes mostly from the inclusion of the N — and 
N = 1 states (first band of negative-parity excited states) , indicating that the ground 
state nucleon couples more strongly to states in these bands. Here the N = line can 
be thought of as a starting point for this work, since the best of the previous studies 
of this state include only the intermediate states included in this band, or restricted 
the intermediate mesons to only the pion. The addition of the N — 2 and N — 3 
bands of states changes the bare mass by a very small amount indicating that the 
sum over intermediate states has converged and a stable solution has been reached 
for the ground state nucleon. This is a very important result, and it will be shown 
that a lack of convergence can greatly affect the final results of such calculations. As 
will be seen below, the impact of the different bands of states varies with the initial 
state studied, and the inclusion of the N = 2 and = 3 band baryons is important 
for other states. In the case of the nucleon, these bands were added for consistency. 

An additional item that needs explanation is the presence of multiple solutions 
for Eq. 5.1 when all intermediate baryon states up to the N — 3 harmonic oscillator 
band are included (solid line in Figure 5.1). This stems from the presence of decay 
thresholds and their effect on the size of the self energies. The locations of some 
ground state thresholds are labeled on the figure, but others cannot be identified 
due to the large number of B'M states included. It is possible that these threshold 
effects could be 'dampened' by the inclusion of the widths of intermediate particles 
as an imaginary part in the energy denominator of Eq. 2.5, but that is a higher-order 
effect which remains to be investigated. In selecting a favored solution for the nucleon, 
studies of the dependence of its self energy on the wave functions have shown that the 
second and third solutions sometimes vanish, but the first solution is always present. 
Therefore it seems prudent to retain only the lower bare mass for the nucleon until 
more consistent results are obtained for the other masses. 
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Self-energy contributions to the mass of the Nucleon 
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Figure 5.1. Sum of the bare energy and self energies as a function of the bare energy 
for the Nucleon ground state with — 0.55 and a — 0.5 GeV. 



Figure 5.2 and Table 5.3(b) show similar results for the ground state A. Again 
the evolution of the bare masses is clear and once more the N — and N — 1 band 
baryon states are seen to make the largest contributions to the self energies. Here the 
convergence is even more apparent, as the addition of the N — 2 and N — 3 band 
states only resulted in an overall downward shift of the E + 11b{E) curves with very 
little movement along the energy axis. The almost vertical slope of the last three 
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lines clearly indicates that not only has the sum over intermediate states converged, 
but that the inclusion of additional states would have an insignificant impact on the 
final bare mass. 

Note that the importance of the labeled thresholds in the case of the A is different 
than those for the nucleon, revealing some of the differences in the internal structure 
of these two states. This refiects the results of similar strong decay calculations which 
predict where to look experimentally for resonances by highlighting strong coupling 
to certain decay channels over others. These calculations also explain why some 
resonances remain 'missing'; they couple weakly to experimentally accessible decay 
channels. 

Combination of the results of Figures 5.1 and 5.2 [also see Table 5.3(b)] shows 
that if only intermediate ground state baryons (A^ = band) are included, the A'"-A 
splitting is roughly 290 MeV. When states in the = 1 band arc included the 
splitting is reduced to about 140 MeV, and that result remains mostly unchanged 
by the addition of baryon intermediate states in the N — 2 and N — 3 bands. 
This agrees well with the expectation from other models (see Ref. [27] where the 
authors find that within their model, 2/3 of the A^-A mass splitting comes from 
one-gluon exchange effects, with the remaining third coming from pion-exchange) 
that a substantial portion of the A^- A splitting should come from a source other than 
the quark-quark residual interactions, in this case a difference in self energies due to 
all B'M intermediate states. This result will be shown to hold despite changes to 
the wave functions from variations in the quark residual interactions. 

5.1.2 Non-Strange L = 1 Negative-Parity Baryons 

The results for the L = 1 negative-parity states are presented in Figures 5.3 
through 5.6 and Table 5.3(b). States with same quantum numbers are shown together 
to facilitate comparisons. 
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Self-energy contributions to the mass of the Delta 




E (GeV) 

Figure 5.2. Sum of the bare energy and self energies as a function of the bare energy 
for the A ground state with — 0.55 and a — 0.5 GeV. 



Figure 5.3 shows the results for the spin-partner states. It can immediately 
be seen that the bands of intermediate states have a different impact on these states 
(and in fact on all the L — 1 negative-parity states) compared to the situation 
with the Nucleon and A shown above. Here, intermediate baryon states up to the 
N = 2 band make sizeable contributions but A*" = 3 states only change the results 
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Self-energy contributions to tine masses of the 

Ni/2-(1535) andNi/2-(1650) 




E (GeV) 

Figure 5.3. Sum of the bare energy and self energies as a function of the bare energy 
for a) A^|"(1535) and b) A^|~(1650) with a, = 0.55 and a = 0.5 GeV. 
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marginally. This indicates that the results have converged and that all intermediate 
baryons up to and including the N — 2 band states are required for convergence. 

If only ground-state baryons are included, the splitting between the two 
states can be seen to be roughly 20 MeV, while it can be observed to grow to roughly 
100 MeV with the inclusion of the A'^ 1 band baryons. Further addition of the 
N — 2 and N — 3 band states brings convergence and closes this gap to about 5 
MeV. This illustrates the wide difference in results that can be obtained if the set of 
intermediate states is not large enough to attain convergence. 

Figure 5.4 shows similar results for the pair. Again the splitting between 

these states induced by these self-energy effects varies from 25 MeV to 170 MeV 
depending on number of intermediate states included. Note that in the case of the 
A"! (1520), although the N — 3 band intermediate states were not required for 
convergence, addition of these states resolved the multiple solution problem. This 
validates the considerable extra effort required to include such a large number of 
intermediate states. Note again the differences in the threshold pattern between the 
two states hinting at how differently these states couple to the various intermediate 
states. 

Figures 5.5 and 5.6 complete the set of non-strange L = 1 negative-parity baryon 
states with the Ai"(1620), A|"(1700), and Ar|"(1675) states. The two A states 
could actually be seen to be spht by a fair amount if only the N — and A" = 1 band 
intermediate states are included, in a way mimicking spin-orbit splitting. However, 
with the addition of the other two bands of states, the A(1620), which was up to 
that point much lighter than its spin partner, becomes almost degenerate with the 
A (1700) and even heavier by roughly 12 MeV. Of course these are states that are 
known to be affected strongly by spin-orbit interactions between the quarks (which 
are not included in this work), so the ordering is, at this time, inconclusive. This 
stresses once again the impact of the choice of intermediate states on the final results. 
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Self-energy contributions to tine masses of the 

N3/2-(1520) and N3/2-(1700) 
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Figure 5.4. Sum of the bare energy and self energies as a function of the bare energy 
for a) Ar|"(1520) and b) A^|~(1700) with a, = 0.55 and a = 0.5 GeV. 
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Self-energy contributions to tine masses of tlie 

Ai/2-(1620)andA3/2-(1700) 




E (GeV) 

Figure 5.5. Sum of the bare energy and self energies as a function of the bare energy 
for a) A|~(1620) and b) A|~(1700) with a, = 0.55 and a = 0.5 GeV. 
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Since the graphs show robust results, we are confident that, for this model, the final 
ordering is the correct one prior to the inclusion of spin-orbit effects. 

5.1.3 = 2 Band States 

Finally, to illustrate a case where convergence is not yet achieved. Figure 5.7 is 
included to show the results for the Roper resonance A^l (1440), with these same 
set of intermediate states. Not only do we not have a unique solution, but it is 
clear from the difference between the bare mass for N < 2 and N < 3 baryon 
intermediate states included in the sum that the addition of more intermediate states 
could significantly change the final result. This indicates the need to extend the 
summation over intermediate states to include at least = 4 band positive parity 
baryon states. Note that the effect of the A'" = 1 band intermediate states is minimal, 
in contrast to the different situations shown earlier. As expected, this indicates that 
as states from higher harmonic oscillator bands are studied, inclusion of a large set 
of intermediate states will be required, with a decreased impact of the lower band 
states and an increased impact of the higher band states. 

5.2 Hamiltonian vs. Self Energies 

As mentioned previously, the splittings between the states resulting from the 
differences in self energies are expected to be comparable in size to those that 
arise from residual interactions between quarks. As a consequence, a self-consistent 
calculation requires that those interactions be adjusted, and with them the wave 
functions, to account for the additional splittings. A priori, we do not know how 
to modify the Hamiltonian to get the desired result of agreement between the bare 
energies and the modified qqq spectrum, as each term affects both the model masses 
and the baryon wave functions. The latter affect the corresponding strong decay 
matrix elements and hence the size of the self energies. To better understand this 
process, the splittings of the states under consideration are examined using several 
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different values of the parameters listed in Table 3.1. Since the best value of the 
^Pq decay strength parameter 7 is affected by changes in the wave functions, it was 
refitted each time to reproduce the strength of the A Nn decay calculated with 
the model of Ref [15]. Additionally, the harmonic oscillator parameter a is chosen 
on a coarse grid to be such that the masses are roughly minimized, with the ground 
state A being at or near its physical mass. 

Four different cases are presented below; first, all quark-quark residual inter- 
actions are turned off leaving only the confining potential to act between quarks; 
second, the contact part of the one-gluon exchange hyperfine interaction is included 
but at about half the strength of the value used in Ref. [15]. Note that although the 
value of as is only marginally lower (see Table 3.1), lowering the value of (Tq has the 
effect of increasing the size of the quarks thereby reducing the strength of short-range 
interactions. In the last two cases, the tensor part is also included (proportionally 
to the contact interaction) and results are presented for two different values of the 
oscillator parameter a. Each table shows the experimental masses of the states, 
the model mass obtained from diagonalisation of the Hamiltonian, and the bare 
masses, extracted from graphs similar to those presented above, solving Eq. 5.1 for 
progressively larger sums of intermediate states. The presence of multiple solutions 
is indicated by a range of energies, or an additional value in parentheses when one 
answer was favored. 

When wave functions with no residual quark-quark interactions are used to 
calculate the self energies, this results in the range of bare masses shown in Table 5.1. 
As can be expected, model states are mostly degenerate but this still results in a 
splitting between the ground state N and A of about 150 MeV due to the difference 
in self energies from the B'M loops. This sphtting comes from flavor and spin 
structure differences between those states, and from differences in how each state 
couples to the various intermediate states included in the sum. It was verified that if 
the masses of all baryons and mesons are assumed to be degenerate, and only ground 
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Table 5.1. Bare masses (GeV) for = 0.0 and a = 0.4 GeV, no residual quark 
interactions. 



State 


Model 


N = 


< 1 


< 2 


N <3 


(Expt. 


Mass) 


Mass 










1 4-, 

[7Vi+] 


(0.938) 


1.230 


1.858 


2.358 


2.367 


2.387-2.492 


[Af] 


1.232) 


1.232 


2.132 


2.508 


2.525 


2.538 




i.ooo J 


1.545 


1.762 


2.500 


2.767 


2.783 


mi 


(1.650) 


1.546 


1.800 


2.487 


2.767 


2.783 


mi 


;i.620) 


1.546 


1.812 


2.467 


2.767 


2.787 


mi 


(1.520) 


1.546 


1.758 


2.400 


2.787 


2.800 


[iVf] 


(1.700) 


1.547 


1.850 


2.537 


2.767 


2.783 


[Af] 


1.700) 


1.546 


2.037 


2.550 


2.878 


2.800 


[All 


(1.675) 


1.547 


2.042 


2.525 


2.800 


2.817 



state intermediate baryons are included in the sum, the splitting between the N 
and A ground states disappears, thereby verifying in this model Zenczykowski [6]'s 
statement about the minimum number of intermediate baryon and meson states to 
be included to reach this symmetry limit. 

Table 5.2 shows the changes created by the inclusion of the contact part of the 
hyperfine interaction (with no tensor or spin-orbit interactions). As mentioned above, 
the contact interaction is roughly half the strength of that used in Ref. [15] and 
the following papers. In the model, the degeneracy is now lifted between the spin 
partners A'"-A, the N^" states, and states. Interestingly, the spin-orbit partners 
A| and A| states remain degenerate. Note that spin-orbit interactions are not 
included in their wave functions. It is important to note that after the addition of all 
intermediate B'M states, the A^-A sphtting remains roughly 150 MeV, demonstrating 
a stable result for these states. Note also that the ordering of the states is seen 
to change as more intermediate states are added, emphasizing the point that not 
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Table 5.2. Bare masses (GeV) for = 0.55, a = 0.4 GeV, with hyperfine contact 
interactions only. 



State 


Model 


N = 


N < 1 


N <2 


N <3 


(Expt 


. Mass) 


Mass 

J.V j.c*oo 










[iVf 


(0.938) 


1.081 


1.812 


2.342 


2.358 


2.375 


[Af] 


(1.232) 


1.232 


2.112 


2.500 


2.508 


2.533 


[Nf 


(1.535) 


1.505 


1.742 


2.108 


2.775 


2.787 


m~ 


(1.650) 


1.588 


1.787 


2.475 


2.750 


2.758 




(1.620) 


1.568 


1.787 


2.200(2.425) 


2.775 


2.787 


[ivf- 


(1.520) 


1.505 


1.750 


2.0625 


2.612 


2.700-2.787 


[ivf- 


(1.700) 


1.588 


1.787 


2.358(2.512) 


2.758 


2.775 


[Af-] 


(1.700) 


1.568 


2.033 


2.358(2.525) 


2.787 


2.792 


[Nf 


(1.675) 


1.588 


2.037 


2.100(2.500) 


2.733(2.775) 


2.787 



all relevant effects have been included by previous calculations which include only 
N — intermediate states. 

Tables 5.3 show the results (previously illustrated in Figures 5.1 through 5.6) 
when both the contact and tensor parts of the hyperfine interaction are included 
in the Hamiltonian. Here the tensor interaction has the strength required from a 
consistent nonrelativistic limit of one-gluon exchange. Each table reflects a different 
value of the harmonic oscillator parameter a (0.4 and 0.5 GeV). Again the A^-A 
splitting is unaffected by the changes, as the wave functions for these two states 
arc essentially unchanged by this change in the basis states. The splittings in the 
bare masses of other states are not strongly affected by the change in a. The model 
masses are minimized with a value of a —0.5 GeV, so the results using this basis are 
preferred. 
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Table 5.3. Bare masses (GeV) for ag — 0.55 with hyperfine contact and tensor 
interactions. 



a) a = 0.4 GeV. 



State Model 
(Expt. Mass) Mass 


iV = 


N < 1 


N <2 


N <3 


[Ar|+](0.938) 1.082 


1.812 


(2.175)2.342 


2.358 


2.375 


[A| J(1.232j 1.232 


Olio 

2.112 


(2.358)2.500 


2.508 


o o c 

2.525 


[iVi"](1.535) 1.500 


1.745 


2.105 


2.712 


2.735 


[N^ J(1.650j 1.572 


1.783 


2.412 


2.737 


2.758 


[Ai~](1.620) 1.570 


1.783 


2.200(2.412) 2.725(2.775) 


2.787 


[iV|"l(1.520) 1.506 


1.750 


2.062 


2.600 


2.650-2.785 


[iV|"l(1.700) 1.606 


1.787 


2.350(2.512) 


2.775 


2.812 


fA|"l(1.700) 1.569 


2.375 


2.350(2.530) 


9 787 


2.795 


[A^r](1.675) 1.584 


2.037 


2.100(2.492) 


2.787 


2.800 




b) 


a = 0.5 GeV. 






State Model 
(Expt. Mass) Mass 


N = {) 


TV < 1 


N <2 


N <3 


[iVi^](0.938) 1.082 


1.850 


2.367 


2.375 


2.392(2.500) 


[A|+](1.232) 1.232 


2.137 


2.508 


2.517 


2.537 


[Ari"](1.535) 1.500 


1.762 


2.375 


2.737 


2.758 


[Ar|"](1.650) 1.572 


1.783 


2.475 


2.742 


2.762 


[A| ] (1.620) 1.570 1.800(2.025) 2.467 


2.800 


2.812 


[iV|"](1.520) 1.506 


1.762 


2.082(2.367) 


2.637(2.793) 


2.800 


[iV|"](1.700) 1.606 


1.787 


2.537 


2.812 


2.825 


[A|"] (1.700) 1.569 


2.050 


2.558 


2.787 


2.800 


[A^|~](1.675) 1.584 


2.050 


(2.100)2.512 


2.787 


2.800 
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5.3 The Spectrum 

Finally this section is brought to a close with the presentation of three figures 
showing the progression of the relationship between the model masses and the bare 
masses required to reproduce the masses of the states extracted from data analyses. 
In each figure the two different mass scales have been adjusted so that the bare and 
model masses of the ground state A coincide. The bare masses shown include all 
intermediates B'M states for each set of parameters. 

When all residual interactions between quarks have been removed, the spectrum 
of the states studied appears as shown in Figure 5.8. As stated above and seen here, 
model masses are degenerate by oscillator band. The inclusion of self energy loops 
requires different bare energies to fit the physical masses of the N and A ground 
states, and reduces the splitting in the bare energies between oscillator bands. The 
Hamiltonian therefore produces model masses in bands that are too far apart at this 
point. 

As seen in Figure 5.9, the addition of reduced-strength contact interactions 
between quarks induces configuration mixing in the wave functions and lifts the 
degeneracy between spin partners. States split by other type of interactions remain 
mostly unchanged by this addition. The effect of the self energies on the bare 
masses extracted from the physical masses also changes, although the A^-A bare 
mass sphtting is stable at roughly 150 MeV. The order of the two states, on 
the other hand, is reversed with the bare mass of the predominantly spin-1/2 state 
heavier than that of the predominantly spin-3/2 state, and the splitting is larger. 
The splitting between the bare masses of A^| states appears to be reduced by the 
inclusion of loops but the presence of multiple solutions blurs the picture somewhat. 

The addition of the tensor part of the hyperfine interaction changes the spectra 
once again, as shown in Figure 5.10. The modifications to the Hamiltonian close 
slightly the splitting between the model masses of the states, and open it for 
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the A^l states. On the other hand the required bare energies of the states are 
almost degenerate, and those of the states are separated by about 25 MeV. It 
is expected that the mixing between states with same quantum numbers due to self 
energy loops will widen the gaps in both cases, so this picture is not yet complete. 
The splitting between oscillator bands is still larger for the model spectrum than for 
the bare mass spectrum but the overall agreement between the spectra remains fairly 
good. 

In closing, it is important to mention the extensive computational work required 
to obtain the results presented in this Chapter and which comprised the most time 
consuming part of this work. The code used to produce the analytical form of 
the momentum dependence of the strong decay vertices was entirely done using 
the symbolic manipulator Maple. The code is based on the general method of 
Roberts and Silvcstrc-Brac [26] and was thoroughly tested by reproducing the large 
number of decay amplitudes found in several published papers by Capstick and 
Roberts [15, 22-23]. The analytical expressions produced by the Maple code were 
subsequently translated into the programming language C, and then included in 
the code which numerically calculates the principal part of the loop integration 
using algorithms such as Gauss-Laguerre and Gauss-Legendre quadratures. These 
extensive calculations would not have been possible without access to the FSU 
Physics Department Computing Cluster. The computational time currently required 
to obtain the results listed in just one of the tables presented above is of the order 
of 5 days of full-time computing for an average of 10 nodes in the cluster. This does 
not, however, reflect the many months of intensive computing work done prior to 
this 'step' when several of the analytical components of the decay amplitudes were 
computed and stored. The interested reader is referred to Appendix C for more 
details about the computational methods used during this project. 
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Self-energy contributions to the mass of the N5/2-(1675) 




Figure 5.6. Sum of the bare energy and self energies as a function of the bare energy 
for A^|"(1675) with a, = 0.55 and a = 0.5 GeV. 
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Self-energy contributions to the mass of the Roper 




Figure 5.7. Sum of the bare energy and self energies as a function of the bare energy 
for A^|^(1440) with a, = 0.55 and a = 0.5 GeV. 
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Figure 5.8. Comparison of the spectra of tlie bare masses required to fit the physical 
masses of the states shown and model masses obtained from a Hamiltonian with no 
residual quark-quark interactions. Here as = 0.0 and a = 0.4 GeV. 
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Figure 5.9. Comparison of the spectra of tlie bare masses required to fit the physical 
masses of tlie states shown and model masses obtained from a Hamiltonian with 
hyperfine contact interactions only. Here = 0.55 and a = 0.4 GeV. 
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Figure 5.10. Comparison of tlie spectra of tlie bare masses required to fit the 
physical masses of the states shown and model masses obtained from a Hamiltonian 
with hyperfine contact and tensor interactions. Here as = 0.55 and a = 0.5 GeV. 
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CHAPTER 6 
CONCLUSION 

Within the context of the relativized quark model, a comprehensive study has 
been carried out of the effects of baryon-meson intermediate states, in the form of 
self-energy loop corrections to the mass of several baryons. The baryon states whose 
self energies have been studied include the nucleon and Delta ground states and 
the first band of negative-parity excited baryons. The ^Pq decay model is used to 
obtain the analytical form of the momentum dependence of the baryon-baryon-meson 
vertices needed in the loop calculations. The model is modified to take into account 
the size of the constituent-quark pair-creation vertex. Wave functions generated from 
a Hamiltonian including reduced-strength one-gluon-exchange interactions between 
quarks are used to calculate the self energies. Since masses play a crucial role in the 
size of the self energies, physical masses are used where known, and model masses [15] 
used otherwise, for the intermediate baryons and mesons. The bare energy of the 
initial baryon is determined self-consistently by calculating the self energies for a 
range of bare energies, then finding the solution to Eq. 5.1 when Mb = Mphysicai for 
each initial baryon studied. 

As demonstrated by Zenczykowski [6] , a minimum set of baryon-meson interme- 
diate states is required to recover the SU(3)/xSU(2)spin symmetry hmit, while Brack 
and Bhaduri [9] showed that intermediate baryon states up to at least the second 
band {N — 3) of negative-parity excited states must be included in order for the sum 
over intermediate baryon-meson states to converge. It is clear that truncation at only 
ground state baryons or only pseudo-scalar mesons leads to physically meaningless 
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results. The present work therefore uses for the first time a complete set of spin-flavor 
symmetry related baryon-meson intermediate states, while at the same time including 
excited baryon states up to the N = 3 band, to insure the convergence of the sum 
of intermediate states. It is shown that drastically different answers are obtained 
for some states if this sum is indeed truncated at smaller sets of intermediate states, 
and that convergence is reached for all the external baryon states considered in this 
work. It is also shown that a larger number of intermediate states is required to reach 
convergence as the initial state becomes more highly excited. For example, the set of 
intermediate states included in this work is found to be insufficient for initial states 
in the first {N — 2) positive-parity excited state band, such as the Roper resonance. 
This and other states in the same harmonic oscillator band will require the inclusion 
of baryon intermediate states up to at least the second {N — 4) positive-parity band. 

The existence of decay thresholds and their effects on the self energies are shown. 
Since physical masses are used where available, actual thresholds can be identified 
on the graphical version of some of the results. Their presence creates an oscillatory 
pattern in the curves for the sum of the bare energy and the self energies, that may 
lead to multiple solutions to Eq. 5.1. Such situations become less frequent as the 
sum over intermediate states is expanded and convergence is reached. 

In this model, it is found that roughly half of the splitting between the nucleon 
and Delta ground states arises from self energy loop effects, the other half coming 
from residual quark-quark interactions. Changes in these interactions have very 
little impact on this result since they do not affect the wave functions of these states 
very strongly. The effects of the same set of intermediate states on the spectrum 
of L = 1 negative-parity excited states is also examined, and it is found that the 
resulting splittings are sensitive to configuration mixing in the baryon wavefunctions 
caused by residual interactions between the quarks. Additionally, some of these 
states are expected to mix further due to off-diagonal terms in their self energies. 
Fairly good agreement is found between the spectrum of bare masses produced by 
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the inclusion of a large set of baryon-meson intermediate states and the spectrum of 
model masses obtained from a Hamiltonian with hyperfine and contact interactions 
between quarks. The overall shift between the two spectra of negative-parity states 
could be attributed to a problem with the string tension, which has been shown [33] 
to be renormalized by the presence of self-energy loops. Spin-orbit interactions are 
also expected to play an important role in changing both model masses and bare 
masses for the negative-parity excited states therefore the results of this work can 
be considered a significant step toward an understanding of these states but work 
remains to be done. 

It is important to note that the model used in this work not only self-consistently 
incorporates the effects on the properties of the quarks of self-energy loops, but also 
the effects of one-boson-exchange mechanism on baryon properties, since the created 
anti-quark can merge with any quark from the initial and final baryon to produce the 
intermediate meson. This point is illustrated in Figure 6.1. Calculations which treat 
only short distance interactions between quarks from the high momentum transfer 
limit of the last two diagrams on the right hand side of Figure 6.1 can be expected 
to neglect important effects from the self-consistent evaluation of quark self energies 
and exchange effects due to all mesons. They also do not properly take into account 
the two-hadron nature of the intermediate state at low momentum scales. 

As mentioned before, and as it should also be apparent from the results of this 
work, conclusions made in prior works about spin-orbit forces in negative-parity 
excited baryon states were premature. Any calculation not including a complete 
set of spin-flavor symmetry related baryon-meson intermediate states, and excited 
baryon states up to at least the N — 3 band, cannot claim to have complete results. 

6.1 Outlook 

After extensive work on this project there still remain many unanswered questions 
and unaddressed concerns. One important asset of the computational tools assembled 
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Figure 6.1. Self-energy model includes both loop effects and one-boson-exchange 
effects in a self-consistent calculation 



for this project is their ability to be extended to larger sets of initial and intermediate 
hadrons given sufficient time and computing resources. Below, a few of the projects 
lined up behind this one are touched upon. 

Since the framework is now in place, the next step is extension to the strange 
sector. The ground state A and E baryons as well as some of the experimentally 
better known L = 1 negative-parity band A and S states are already under study. 
Experimental data for several of these states has recently improved, but agreement 
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between the analyses of this experimental data and model predictions for resonance 
parameters is far from perfect. It will be interesting to find out how the baryon- meson 
loop induced self energies affect the sphttings within this set of states. 

As pointed out earlier, the inclusion of at least N = A band baryon intermediate 
states will be required to bring about convergence of the intermediate state sum 
ior N — 2 band initial states. This should shed some light on the nature of some 
controversial and hard to model positive-parity states such as the Roper resonance. 
The computational tools used in this work are easily extended to such a calculation. 

Mixing due to self energies must contribute to the splitting between baryons with 
same quantum numbers. Therefore studies of the mixing between the Nucleon and 
the Roper, the A^|"(1535) and A^|~(1650) states, and the A^|~(1520) and A^|~(1700) 
states will be required to further understand the impact of self energy loops on the 
splittings of these pairs of states. 

The addition of spin-orbit interactions needs to be further explored. A prelim- 
inary study shows that the addition of spin-orbit interactions in the interactions 
leading to the wave functions used for the decay vertices can change the splitting 
between some states and even reverse the ordering of some of the states. More work 
is clearly needed in this area before conclusions can be drawn, as both the mixings 
due to self energies and spin-obit interactions strongly affect the masses of these 
states. 

Even with these cautions, the present calculation goes far beyond anything 
previously available for the negative-parity non-strange baryons and shows, for 
the first time, results with a set of intermediate states large enough to achieve 
convergence. It also demonstrates for the first time the sensitivity of the self energies 
to the mixings caused by various components of the quark residual interactions. 
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APPENDIX A 



THE WAVE FUNCTIONS 

There are four components in each hadron wave function: color (C), flavor (0), 
spin (x), and spatial (-0) wave functions. 

A.l Color 

The meson color wave function is found by the rules of SU (3)c for direct products 
of quarks {q) and antiquarks {q) carrying color charges or3(8)3 = l®8 and must 
be singlet to be an observable therefore it is 

= E ^ q^ Qi- (A.l) 

Similarly, the baryon color wave function is found by the same rules but for three 
quarks or3(8)3(8)3 = l08®8'®lOtobe the totally antisymmetric singlet 
combination 

C|= E -^e^jkqUiql (A.2) 
i,j,k=i V6 

A.2 Flavor 

The flavor wave functions for the baryons and mesons included in this calcula- 
tion are shown in Table A.l. They are obtained from the irreducible representations 

of SU (3) F giving a flavor nonet each for the pseudoscalar mesons (tt, K, K, t], t}') and 
the vector mesons (p, K* , K*,uj, (py. 

For baryons, we follow the convention used in ref. [15] and adopt a generalized 

uds basis that only symmetrizes the product (pxi^ in identical quarks. This removes 

^Note that since the (p mesons couple weakly to non-strange baryon states (such decays are OZI 
suppressed), they are not included. 
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Table A.l. The baryon and meson flavor wave functions as a function of their 
isospin projection (J^). 



State +3/2 +1 1/2 -1/2 -1 -3/2 

N uud ddu 

A uuu uud ddu ddd 

A -^{ud — du)s 

S uus -^{ud + du)s dds 

TT —ud — dd) +du 

p —ud 71 ~ 

X —MS +SU 

-ds -sd 

uj -j^{uu-\-dd) 



the need for symmetrization between the u and d quarks in the spatial wave function 
making it manageable to work with states with up to ITiuj in the harmonic oscillator 
spectrum. Note that the baryon flavor wave functions are all either symmetric or 
antisymmetric under the interchange of quarks one and two. 

A. 3 Spin 

The total spin of two spin-i particles can be either zero or one giving us the 
pseudoscalar and vector mesons respectively. 

The total spin of the three spin-| particles can be either | or | so that as a 
complete set of spin wave functions x we can choose 

xfs = im) , etc. (A.3) 

2 2 

xff = 4(inT)-|m)) , etc. (A.4) 

2 2 ^J2 
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Figure A.l. Relative coordinates p and A. 

xtt = -4(im) + im)-2im)) , etc. (a.s) 

2 2 1/6 

(We show only the top state of a JM multiplet; other wave functions follow the 
Condon-Shortley convention). Note that the baryon spin wave functions are also 
either symmetric or antisymmetric under interchange of the first two quarks. 

A. 4 Space 

Finally, for the spatial wave functions ^ we take functions with definite total 
L = Ip + 1a made from a Clebsch-Gordan sum of harmonic oscillator wave functions 
in the two relative coordinates 

P = 75(^1-^2) (A.6) 

and 

A ^ -J=(ri + r2-2r3) (A.7) 
of the three body problem (see Figure A.l). These are 
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LMriplpTi^l^ 



(A.8) 



where the Ln ^ (x) are the associated Laguerre polynomials 



1\ ^2m 



m=0 



n — m m\ 



(A.9) 



(half-integral factorials are defined by the F function), and the normalization 
coefficient Mni is defined by 



nl 



2n\ 



(A.IO) 



Putting all the elements together, the wave function is then expanded in a set of 
states of the form 



) = Ca^ C{L S Ml J -Ml] J M)^ 



Ml 



(A.ll) 



The entire wave function is now explicitly antisymmetric under the exchange of 
quarks one and two. 
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APPENDIX B 



TRANSITION AMPLITUDE 



The final form of tlie transtition amplitude is 
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Here 



with 



and 



{ Sa Lp -Sa 1 f Lp Sp Jp \ [ Lp *6 1 r Lp Sp Jp 1 

14 Ja La\\\ Sa Sa } \ h Jb L, j \ ^ Sj 



Sb Lp Sb 

2 "a i^a J y Jb Lb 

[-iY+ta+Jc-u-Scj:(^j^BC)n{ABC) 



xE(-i) 

She 



Jp 1/2 Sb 
1/2 1/2 S, 

Sa 1 She 



Lbc 



X 



Ef 2 / 'Sa 4 Ja \ ( 
. I L Sbc 1/1 



^bc J^bc Jbc 
e Ja L 



Sb h Jb 

Sc. Lc Jc 
Sbc Lbc Jbc 

e{eb,L„Lbc,e,ea,L,ko). (b.i) 



Ja — La + Sa — 4 ~l~ ^ai 

La = + Lp„ = 4 + Lp^, 

Sa = S,„ + l/2, 



Sa — Jpa 



1/2 = 



1/2, 



(B.2) 



(B.3) 



(B.4) 



with similar definitions for B. The first four 6 — j symbols of Eq. (B.I) are necessary 
for transforming from the usual angular momentum basis for the baryons, given by 
Eq. (B.3), to the basis of Eq. (B.4), which is the more convenient one for evaluating 
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the transition amplitude. L, L^c and S^c are internal summation variables, and 
!F{ABC) is the flavor overlap for the decay. 

The purely "spatial" part of the transition amplitude is 



X E 

^12,^5/6,^7/ 



L 





^'1 


4 " 




"4 


<3 


1 " 










4 




4 


12 


4 


Lbc_ 




4 


4 


L 



X 



^6 -f' -f/bc 



D^12 Ejfe D^6 D^7 D^8 

^hh ^Ui2 ^t'^e^ ^iaU ^e'^e'^ 



X j^^l +^2 +^3 +€4 +2 A+i/y Q2ij.+u-ei -ti-tz-U 



^;+£'2 + 4 + 4 + 2At + i/ + r 



(B.5) 



In this expression, is a normalization coefficient that results from writing a 
single component of the wave function of A as 

*LMn,^,nx€A (Pi, P2, Ps) = T]{AJ^fl'^ ^ < 44mM - m\LM > 



X 



For proper exchange symmetry among the quarks. A' — -^A^ and 



1 1 
Pp = 2 ~ P2) ' = 3 + P2 - 2P3) ■ 



(B.6) 



(B.7) 



77 is a phase factor that arises from calculating the Fourier transform of the 
configuration space wave functions, and has the value 



(B.8) 



With these definitions, Na = A^^+'^/^A/'„a«a' "with J\fni previously defined in A. 10, 
^£+1/2 ^ while the Yg^ are the usual spherical harmonics. 
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i7 is a Jacobian factor needed to convert from the basis used in evaluating the 
space factor e in Ref . [26] , to the basis used in the evaluation of the wave functions 
used for exphcit calculation of the decay amphtudes. The wave functions of Ref. [21] 
use 

Pp = ^ (Pi - P2) ' Pa = ^ (Pi + P2 - 2P3) , (B.9) 

so that both the Jacobian factor mentioned above, as well as a redefinition of the 
gaussian parameters of the wave functions, are required in order to use the wave 
functions of Ref. [21] with the above expression for the decay amplitude. 

The factor TZ of Eq. (B.l) is obtained as the overlap of the wave functions in the 
p coordinates in the initial and final baryon. Since in the model used here quarks 
1 and 2 are spectators {ip^ = ip^, Sp^ = Sp^^, Jp^ = JpJ, and the basis is fully 



orthogonalized {a is the same in the initial and final baryons, so that n^^ 



n 



pbJi 



this 



overlap is always unity. In addition, this means that the Jacobian discussed above is 
only necessary for the transformation in p;^. 

The T,x,ij.,„ Dx^,,{oJi,u;2,x)It,(£5,£6,£7,£8; L) term arises from writing (here Qa = 
Pa„, with a similar definition for the daughter baryon) 



(B.IO) 



When the substitutions = + q, q^ = (a; — uji)k + q, qc = {x — 002)^ + q are 
made, and the integrals over k and q are evaluated, the expression above results. 
The full form of the Dx^u does not provide additional information so it omitted. 
In Eqs. (B.l) and (B.5), 



(B.ll) 
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b 


c 
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b 
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= cfghi < 
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e 




.9 
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i 




.9 


h 





where < 



a b 
d e 
, 9 h 



/ > is the 9 — j symbol, and J = y/2J + 1. 
i J 
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In Eq. (B.5) 



1 



+ B^x- coif + C^x- co^f + f{x - If 



2 



(B.12) 



uji and LJ2 are ratios of various linear combinations of quark masses. In general, 

m3 



mi + m2 

UJi = ■ ■ , LJ2 



mi + m2 + mi m^ + m^ 

where the subscripts refer to the quark labels shown in Figure A.l. In addition, 



(B.13) 



C 



47r(2£+ 1)! 



B. 



\ (2£i + l)![2(£-£i) + l]!' 





(B.14) 



and l\ — Li, — £i, £'2 = Lc — £2, '^3 — ^ — '^3: '^4 — La — £4 and the geometric factor 
is 



/2p(4,4,^7,4;^) = (-i)^(2p)!4W8 

^ 4^(4A + l)(p + A)! f2\ 4 e7\f2X 4 4 





''5,(2p + 2A + l)!(p-A)! V 



/2p+i(4, 4, 4, 4; L) = 2(-i)^+i(2p + i)!4444 
^ 4^(4A + 3)(p + A + l)! /2A + 1 4 4 
(2p + 2A + 3)!(p-A)! V 



/2A + 

n 



1 4 4 




4 4 L 
is 4 2A + 



x} 



4 4 i> 
4 -^7 2A 



(B.15) 
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APPENDIX C 



COMPUTATIONAL METHODS 



There are two main computational elements required to accomplish the calcula- 
tion in this work. Recall the main equation of Chapter 2 



where A^bb'm(^) is an expression for the momentum dependent vertex of the strong 
decay B —>■ B'M, and V indicates that only the real part of the integral is evaluated 
by principal-part integration. How these components are computed is described in 
general terms in what follows, along with a sequence of how they are combined to 

produce the results presented in Chapter 5. 



The momentum dependence of the strong decay vertex can be obtained from the 
analytical form of the strong decay matrix element {B'M\T\B) , with T representing 
the ^Po decay operator. The method used to evaluate these matrix elements is 
based on the work of Ref. [26] and Ref. [15], and several of the angular momentum 
techniques can be found in Ref. [28]. 

The symbolic manipulator Maple was used to compute all the components 
involved in calculating each decay matrix element for a large set of decays B — > B'M. 
Subroutines were built to independently calculate each component so that it could be 
individually tested before being integrated into a higher level subroutine, and so that 
it could be totally portable. Basic procedures such as those to analytically calculate 
Clebsch-Gordan, 6 — j, or 9 — j coefficients were implemented, along with generators 




(C.l) 



C.l Strong Decay Vertices 
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of spherical harmonics and associated Laguerre polynomials. Then another series of 
algorithms were devised and coded to handle each component of Eq. 4.11 as defined in 
Appendix B. The very large number of nested summations involved required careful 
and thorough testing as the trouble-shooting became more complex with each 'layer'. 

Due to the symmetry of the wave functions, in this ^Pq model decays involving 
the first or second quark are considered separately from those involving the third 
quark of the decaying baryon. However, the two cases are related by a set of 
coefficients that, in effect, project a set of basis states onto another, or equivalently 
rotate a set of coordinates. The matrix elements can then be calculated in two 
different bases, and then combined after one set of results is transformed. These 
coefficients were adapted for this work from Moshinsky [29-31] brackets, coded, and 
then subsequently mass-produced with the results stored in analytical form. Once 
that was accomplished, it became possible to rapidly analytically compute matrix 
elements for decays from any initial baryon to any baryon-meson state, making this 
code a very versatile and powerful tool. 

Since Maple is a powerful but slow program, and large parts of the analytical decay 
matrix elements are common for specific sets of quantum numbers, matrices of decay 
matrix elements were created between sets of hadron quantum numbers based on the 
expansion of the baryon wave functions. For example, matrices were built for decays 
from — I"*" baryons {e.g. the A baryon) to — |^ baryons {e.g. the nucleon) 
and a pseudo-scalar meson {e.g. the pion). Since the baryon wave functions used in 
this project are expanded to the = 6 level (for positive-parity states and N = 7 for 
negative-parity states) and thus have an average of about 100 components, for this 
example one matrix is 78 x 50, the other 152 x 100, with each element corresponding 
to a strong decay between a pair of basis substates. The process is repeated for 
decays involving vector mesons {e.g. the rho). Given the number of intermediate 
baryon-meson states included in this work, it is an understatement to say that this 
was a computer-intensive endeavor. It took several months to complete, using many 
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nodes on the FSU Computer Cluster. Matrices were finally tested by using them to 
reproduce a large variety of pubhshed decay amphtudes [15, 22-23]. 

The reward for this extensive calculation is the versatihty of the code. Not only 
can it be used in the type of project described in this dissertation but, as mentioned 
above, it can also be applied to obtain any strong decay amplitude within the limits of 
the ^Pq model. Since the decay matrix elements are stored in analytical form, changes 
in the value of any parameter can easily be handled without having to recalculate 
anything, and since wave function expansion coefficients are independent of the decay 
matrix elements, they can also be changed without affecting this part of the code. 

C.2 Numerical Integration 

Since the integrand of Eq. C.l can become extremely complex, it is not practical, 
and often not possible, to try to do the integration analytically, so numerical 
integration schemes were therefore used. Since both real (the initial baryon's energy 
is above the threshold for production of the intermediate baryon-meson pair) and 
virtual (below threshold) decays are encountered, two different routines were used. In 
the case of virtual decays the integrand is always real, so a Gauss-Laguerre quadrature 
routine was used to perform the integration. Above threshold decays imply the 
presence of a pole, therefore a combination of Gauss-Laguerre and Gauss-Legendre 
type quadratures were used to integrate symmetrically around the pole and evaluate 
the principal part of Eq. C.l. 

The integration routines were developed using the C programming language, and 
used modified versions of pre-coded numerical algorithms from Ref. [32]. 

C.3 Overall Scheme 

The numerator of Eq. C.l is obtained via a Maple routine which performs the 
matrix algebra required to combine the matrices of decay matrix elements with 
the appropriate vectors of wave function expansion coefficients. The result is an 
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analytical expression as a function k which is then translated into C code by Maple. 
The result is then used by the integration routine. 

During the matrix algebra part of the calculation, it was found that the amount 
of time required to process analytical matrices ranging in size from 150 MB (for 
decays involving pseudo-scalar mesons) to well over 650 MB (for decays involving 
vector mesons) could limit the number of intermediate baryon-meson states included. 
Since a crucial ingredient of the calculation is the large extent of the sum over 
intermediate states, an additional routine was implemented to 'preprocess' the 
matrices by assigning numerical values to the parameters that were common to a 
set of initial states (based on a given set of baryon wave functions) , which were then 
stored in new matrices for use when needed. This extra step, added at the cost of a 
few days of processing for each new set of matrices, cut the final matrix algebra time 
by a factor of roughly 60%, making it possible to compare multiple sets of results 
within a reasonable time frame. 

For each initial baryon and each baryon-meson combination in the sum over 
intermediate states {^b'm)i the associated integral is evaluated for 200 different 
values of E over the energy range being considered. These values are then tabulated 
and then used to produce graphs such as those presented in Chapter 5. 

Throughout this process, liberal use is made of scripts (both in Maple and Unix) 
to handle the large number of intermediate states and automate the process as much 
as possible. 
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